Strong coupling superconductivity, pseudogap and Mott transition 
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An intricate interplay between superconductivity, pseudogap and Mott transition, either band- 
width driven or doping driven, occurs in materials. Layered organic conductors and cuprates offer 
two prime examples. We provide a unified perspective of this interplay in the two-dimensional 
Hubbard model within cellular dynamical mean-field theory on a 2 x 2 plaquette and using the 
continuous-time quantum Monte Carlo method as impurity solver. Both at half filling and at finite 
doping, the metallic normal state close to the Mott insulator is unstable to d-wave superconductivity. 
Superconductivity can destroy the first-order transition that separates the pseudogap phase from 
the overdoped metal, yet that normal state transition leaves its marks on the dynamic properties 
of the superconducting phase. For example, as a function of doping one finds a rapid change in 
the particle-hole asymmetry of the superconducting density of states. In the doped Mott insulator, 
the dynamical mean-field superconducting transition temperature T^ does not scale with the order 
parameter when there is a normal-state pseudogap. T^ corresponds to the local pair formation 
temperature observed in tunneling experiments and is distinct from the pseudogap temperature. 

PACS numbers: 71.27.+a, 71.10.Fd, 71.10.Hf, 71.30.-fh 



The proximity between a Mott insulator and a su- 
perconductor is one of the most intriguing puzzles in 
condensed matter physics [ ]. Indeed, in a Mott insu- 
lator, strong Coulomb repulsion between electrons is at 
the origin of the phenomenon, while superconductivity is 
usually associated with effective attraction. In half-filled 
band layered organic conductor, pressure induces a first- 
order transition between a d-wave superconductor and 
a Mott insulator. This is a bandwidth-induced transi- 
tion. The maximum superconducting transition temper- 
ature Tc is at the first-order phase boundary [-]. On the 
contrary, in high-temperature superconductors, while su- 
perconductivity emerges upon doping a Mott insulator, 
Tc has a dome shape and disappears before the doping 
driven Mott transition [ ]. In addition, the normal state 
near the Mott insulator exhibits a pseudogap [4]. 

Weak coupling approaches to the simplest model that 
includes screened Coulomb interaction and band struc- 
ture effects, the Hubbard model, show that d-wave su- 
perconductivity can arise as a secondary effect from 
exchange of antiferromagnetic fluctuations [■'3-10]. At 
strong coupling, renormalized mean-field theory [11-13], 
slave particle [14, 1-")] and variational approaches [16, 17] 
also suggest the presence of d-wave superconductivity. 
However, to study both the Mott transition and d-wave 
superconductivity, one must resort to cluster versions of 
dynamical mean-field theory [18-21]. Up to now, results 
have been obtained mostly at zero temperature [22-31]. 
There are also a few results on the transition tempera- 
ture [32-3-'')] but there is no systematic study of the in- 
terplay of superconductivity and pseudogap with both 
bandwidth-driven and doping-driven Mott transitions at 



finite temperature. This is the problem that we solve 
in this Letter by studying the two-dimensional Hubbard 
model with cellular dynamical mean-field theory on a 
plaquette [20, 21] using state of the art continuous-time 
Quantum Monte Carlo method as impurity solver [:^)()- 
:]!)]. Notice that quite generally [40] there is no con- 
tinuous symmetry breaking in two dimensions at finite 
temperature. This is true for d-wave superconductivity 
as well [ ]. However, it is still physically meaningful to 
study the superconducting phase at the dynamical mean- 
field level since the corresponding transition temperature 
indicates where the superconducting fluctuations be- 
gin to develop. Three-dimensional effects eventually al- 
low true long-range order at lower temperature. Compe- 
tition with other long-range ordered phases [42], which 
are infiuenced by many factors including frustration, will 
be considered in future work. 

After we present the model and method, we discuss in 
turn the bandwidth-driven and the doping-driven cases 
before we provide a unified view and discussion of the 
results. 

Model and method.- We consider the two-dimensional 
Hubbard model on a square lattice, 

H ^ -J2t,^clc,^+U J2 (^»T - 1/2) (^^4 - 1/2)-mE 

(1) 

where and Cia create and annihilate an electron of spin 
a on site i, ni„ — c^^Cia, t is the nearest neighbor hop- 
ping amplitude, /i is the chemical potential and U is the 
screened Coulomb repulsion. We solve this model using 
cellular dynamical mean-field theory (CDMFT) [20, 21]. 
This approach takes a cluster of lattice sites, here a 2 x 2 



plaquette, out of the lattice, and embeds it in a self- 
consistent bath of noninteracting electrons. The action 
of the plaquette coupled to the bath reads 



dr 



drVt(^)A(^,/)^(^'), (2) 



where Sc is the action of the cluster and A the hybridiza- 
tion matrix. From now on, the symbol " indicates a ma- 
trix in cluster indices. The hybridization A is determined 
by the self-consistency condition 



A(iw„) 



(3) 



which states that the infinite lattice and plaquette have 
the same self-energy and the same Green's function on 
the plaquette. Here Sc is the cluster self-energy, tc the 
cluster hopping, and G(*^„) = Efc ' 

where k is the superlattice momentum. We solve the 
impurity (plaquette-|-bath) problem of Eq. (2) using the 
continuous-time quantum Monte Carlo method [36, 39], 
which sums all diagrams obtained by the expansion of 
the action of Eq. (2) with respect to the hybridization A. 
For the superconducting state in the cluster momentum 
basis, the cluster Nambu Green's function reads 



Gk{t) = 



(4) 



Gki;{t) Fk{t) 
F+{t) -G_K;(-r) 

where F is the anomalous Green's function. For d-wave 
superconductivity, i^(,r,o) = ^-P'co.tt) is the only nonzero 
component. To determine the parameter space where the 
superconducting phase is allowed by the CDMFT equa- 
tions, we monitor the superconducting order parameter 

*=(i^KO)(T = 0+)). 

Superconductivity and interaction- driven Mott 
transition-First, consider the normal state of the 
half-filled two dimensional Hubbard model. Previous 
work revealed a first-order transition at moderate 
interaction between a correlated metal and a Mott 
insulator [43-45]. As shown in Fig. la, in the {U,T) 
plane there is a hysteresis region (in red or light gray) 
where two mean-field solutions can be obtained. This 
region is bounded by the spinodals Uci{T) and Uc2{T) 
(red lines with triangles) where the double occupation 
shows sudden jumps. The first-order metal-insulator 
transition lies within this region and starts at the critical 
Mott endpoint (C/mit,Tmit) ~ (5.95i, l/12t). 

Next we allow for d-wave symmetry breaking in the 
CDMFT equations and perform scans as a function of U 
for different temperatures. As input seed of the CDMFT 
iterative procedure we use the normal state converged 
solution, and we add a small perturbation in the anoma- 
lous component of the hybridization matrix. We obtain 
a converged superconducting solution, characterized by 
a nonzero close to the Mott transition. No supercon- 
ducting solution is found if we use the metastable insu- 
lating solution as seed. Figure lb shows the order param- 
eter $ for the low temperature T/t = 1/100. Within our 




FIG. 1; (a) Temperature T versus interaction strength U 
phase diagram of the half-filled two-dimensional Hubbard 
model obtained by CDMFT. Three phases can be distin- 
guished: correlated metal, Mott insulator, and superconduc- 
tor. In the normal state, there is a first-order transition at 
finite temperature between a correlated metal and a Mott in- 
sulator, bounded by the spinodals Uc2{T) and Uci{T), defined 
as the loci where the double occupation shows a jump. The 
superconducting phase (blue or dark gray region) is defined 
by the loci where |$| 7^ (filled blue circles) and is dehmited 
by the superconducting transition temperature T^. Extrapo- 
lations to r = are a guide to the eye. On the right vertical 
axis we convert to physical units by using t = 0.35eV. Inset: 
zoom on the superconducting phase, (b) d-wave supercon- 
ducting order parameter <E> as a function of U at half filling 
and for T/t = 1/100. (c ) Density of states p{lo) for U = 5.6t 
and T/t = 1/100 for the normal-state Mott insulator, the 
normal-state metal and the superconductor (dotted, dashed, 
and solid lines, respectively). 



numerical precision, as a function of U, the order param- 
eter exhibits two jumps: one at U{T/t = 1/100) ~ 5.45 
where there is a transition from the metal to the super- 
conductor, and one at Uc2{T/t = 1/100) ss 5.65 where 
the transition is between the superconductor and the in- 
sulator. 

By performing the above procedure for different tem- 
peratures, we obtain the superconducting region in the 
([/, T) plane [blue or dark gray region in Fig. la], defined 
as the region where $ 7^ 0. With decreasing temperature, 
the superconducting phase emerges from the normal state 
metal close to the Mott transition, i.e. for U < Uc2, and 
rapidly disappears below Ud- The largest superconduct- 
ing transition temperature T^{U) occurs, along with the 
largest order parameter, around the first-order boundary 
with the insulator, as in the organics [2]. 

Physically, the CDMFT superconducting transition 
temperature is the temperature below which Cooper 
pairs form within the cluster. Previous work [33] suggests 
that converges to a finite value with cluster sizes up 
to 26 sites. Long-wavelength thermal and quantum fiuc- 
tuations in the magnitude [ ] and phase of the order 
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FIG. 2: (a) Temperature T versus doping 5 phase diagram 
at [/ = 6.2t > C/mit, obtained by CDMFT. Four phases can 
be recognized: in the normal state, there is a first-order tran- 
sition at finite temperature between a pseudogap and a cor- 
related metal, bounded by the spinodals Sci{T) and Sc2{T) 
(up and down triangles respectively). A crossover takes place 
above the critical endpoint {Sp, Tp) and defines the pseudogap 
temperature T* [oi], determined by max dp{Lj — 0)/dT. The 
third phase is the Mott insulator at 5 = (green solid line). 
The fourth phase, the superconducting one, is delimited by 
Tc{5), i.e. the temperature below which |"I>| 7^ 0. Extrap- 
olations to T = are a guide to the eye. Inset: chemical 
potential ^ versus doping 5 = l — na,tT = 1/100 for the nor- 
mal state (triangles) and the superconducting state (circles). 
The jump in the dopings identify the spinodal points between 
the two normal-state metals, i.e. the pseudogap (PG) and 
the correlated metal (CM). The transition is removed by the 
superconducting state: fi{5) does not show any sign of hys- 
teresis, (b) d-wave superconducting order parameter $ as 
a function of doping for temperatures T = 1/64 > Tp and 
1/100 < Tp. On the right vertical axis we convert to physical 
units by using t = 0.35eV. 

parameter [47- .50] will lead to an actual superconduct- 
ing transition temperature smaller than T^. Long- 
wavelength antiferromagnetic fluctuations on the other 
hand can increase T^, as seen in weak-coupling calcula- 
tions [•">-!)]. Competing long-range order would reduce 
or eliminate [ -]■ Nevertheless, informs us on the 
regime of temperature where strong coupling and short- 
range nonlocal correlations lead to pairing. These effects 
lead to a strong d-wave pairing gap in the density of 
states of Fig. Ic. 

Superconductivity and doping-driven Mott transition- 
We turn to the doped Mott insulator. Previously, we 
explored the normal-state phase diagram [52, 53] and 
demonstrated that the first-order transition at half filling 
naturally extends at finite doping, and that it can take 
place between two metallic states: a correlated metal at 
large doping and a pseudogap [51]. Figure 2a shows the 
{5,T) plane at f7 = 6.2t > C/mit- The spinodals 6^i{T) 
and 5c2{T), determined by the jumps in the doping 5 
(see inset), envelop the transition and terminate at the 



critical point (i5p, Tp), which is the extension of the Mott 
critical point away from half filling. The value of {6p,Tp) 
moves to larger dopings and smaller temperatures as U 
increases. At U = 6.2t, Tp is sufficiently large to be acces- 
sible by simulations. Associated with the critical point 
{Sp,Tp) there is a Widom line [54], and the pseudogap 
temperature T*{S) occurs along this line [51]. 

Next, we study the superconducting phase as a func- 
tion of doping. The superconducting order parameter is 
shown in Fig. 2b for different low temperatures. In the 
Mott insulator at zero doping, <I> = and thus there is 
no superconductivity. Upon hole doping, $ increases, 
reaches a maximum for a doping near the normal-state 
first-order transition between the pseudogap and corre- 
lated metal, and, with further doping, decreases. 

By monitoring <&(^) for different temperatures, we can 
construct the superconducting region in the {S, T) plane 
(blue/dark grey region in Fig. 2a). The transition tem- 
perature is higher than the critical temperature Tp, 
and superconductivity eliminates the first-order transi- 
tion of the underlying normal state. Indeed, the (5(/i) 
curve in the inset of Fig. 2a) is continuous. T^ is zero 
at (5 = 0, but it is finite for ^ — >■ 0+ and does not show 
large variation when there is a pseudogap in the under- 
lying normal state. In particular, T^{5) does not ap- 
preciably increase as we approach half filling while the 
pseudogap temperature T* does, showing that the two 
phenomena are distinct, as also found in high- field trans- 
port measurements [55, 56]. With further doping, when 
the superconductivity evolves from a correlated metal, 
T^ decreases and vanishes at large doping. Therefore, 
our results imply that Mott physics causes $ to drop at 
small doping, but does not produce a fall in T^. cor- 
responds to Cooper pair formation within the plaquette. 
We associate T^ to the temperature at which a super- 
conducting gap appears in tunneling experiments [57, 58] 
without long-range phase coherence. Experimentally, in 
the doping range where there is a normal-state pseudo- 
gap, that temperature scale is smaller than T* and larger 
than the actual Tc. The small value of $ suggests that 
the actual Tc of the system will vanish at small doping 
due to competing order [42] or to disorder [59, 60] or to 
long wavelength (classical and quantum) fluctuations of 
the magnitude [4(i] or of the phase [ - ] of the order 
parameter. 

Even though superconductivity eliminates the first- 
order transition in the underlying normal state, signa- 
tures of that transition remain in the dynamics of the 
superconducting state. This is shown by the evolution 
of the density of states with doping in Fig. 3 where the 
solid line is for the superconducting state and the dashed 
line for the normal state. At low doping, superconductiv- 
ity originates from the pseudogap and the superconduct- 
ing density of states inherits its large particle-hole asym- 
metry [34], as found in experiments [(il]; On the other 
side of the transition, at large doping, superconductivity 
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FIG. 3: Low frequency part of the local density of states p{ijj) 
at (7 = 6.2t, T/t = 1/100 for the normal-state and the su- 
perconducting state (red dashed and blue solid lines). For 
S ~ 0.01,0.03,0.04 [panels a, b, c respectively] the super- 
conducting state emerges from the underlying normal-state 
pseudogap metal. It inherits a strong particle-hole asymme- 
try. For 5 ~ 0.06 the superconducting state emerges from a 
correlated normal-state metal, and the density of states, near 
the transition, approximately recovers particle-hole symmetry 
at low frequency. 



emerges from the normal-state correlated metal, and the 
superconducting density of states at low frequency close 
to the normal-state transition is particle-hole symmetric. 
Our contribution is to link the features of the supercon- 
ducting density of states to the underlying normal-state 
first-order transition. 

Superconductivity from Mott physics. -The above anal- 
ysis shows that superconductivity arises by approaching 
the Mott insulator as a function of both the interaction 
strength and the doping. The two routes to create su- 
perconductivity are related, as sketched by the ([/, /i,r) 
phase diagram in Fig. 4. The critical endpoint {fj,p,Tp), 
hidden by the superconducting phase in the (/i, T) plane 
of that figure, is connected to the familiar Mott endpoint 
(C/mit,?mit) at half filling (see dotted line in Fig. 4). 
The latter appears above the superconducting phase. Re- 
cent works at half filling [ ] did not find a direct transi- 
tion between the superconductor and the Mott insulator. 

Previous CDMFT works [25, 27-29, 31, 34, G2] at zero 
temperature reported a doping dependence of the order 
parameter $ similar to the one found here, but the doping 
dependence of could only be surmised. Our contribu- 
tion is to show that does not scale with $(5) when 
a pseudogap is the underlying normal state. remains 
finite as the Mott insulator is approached, implying that 
Mott physics does not suppress even though it sup- 
presses the order parameter. In the region where there 
is a normal-state pseudogap, represents a local pair- 
formation [-57, 58] temperature scale that is distinct from 
both T* and the actual superconducting long-range phase 
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FIG. 4: Schematic temperature - chemical potential - interac- 
tion strength phase diagram based on CDMFT solution of the 
2D Hubbard model. Cut at particle-hole symmetry {/i = 0) 
and at constant U > Umit are shown. Since we set t' — 0, 
the phase diagram is symmetric with respect to fi = plane. 
The first-order transition between a metal and a Mott insu- 
lator in the fi = plane is connected with the first-order 
transition between the pseudogap and a correlated metal in 
the U > C/mit plane [■"■_!, '•'■>]. Tp begins at Tmit. The super- 
conducting temperature T^, delimiting the region where "I> is 
non zero, is also shown. In the phase diagram, the supercon- 
ducting phase emerges from the normal state metal close to 
the Mott insulator. 

coherence Tc- In addition, we find that a classical, not 
quantum, critical point at finite temperature between a 
pseudogap and a correlated metal [51-5.')], continues to 
control the distinct pseudogap physics above T^, even 
though the superconducting phase replaces the normal- 
state first-order transition at low temperature. This find- 
ing has to be contrasted with the quantum critical point 
reported in previous work [()3]. Because those calcula- 
tions were limited to high temperatures, they did not 
detect the normal-state first-order transition. 

The phase diagram as a function of interaction 
strength, doping and temperature that we found shows 
that a transition directly to the superconducting state 
from a Mott insulator is possible at the dynamical mean- 
field level, whether the transition is bandwidth or dop- 
ing driven. Since is finite at infinitesimal doping, the 
transition appears as first-order in both cases. Hence, the 
experimentally observed drop of at low doping must 
come from mechanisms not included here, such as long 
wavelength fluctuations [46-50], competing order [42] or 
disorder [59, 60]. Long- wavelength fluctuations should 
be important near the Mott transition because the order 
parameter decreases rapidly with decreasing doping, con- 
trary to T^. Yet, retains a role as a local pair forma- 
tion temperature [57, 5S] and is distinct from the pseudo- 
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gap temperature T*. For sufficiently large U the super- 
conducting state destroys the underlying first-order tran- 
sition between the pseudogap and the correlated metal, 
but signatures of this transition remain in the dynamical 
properties of the superconductor. 
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